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Abstract. We apply optimal control theory to a tuberculosis model given 
by a system of ordinary differential equations. Optimal control strategies are 
proposed to minimize the cost of interventions. Numerical simulations are 
given using data from Angola. 
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,J^ ' 1. Introduction. Angola has one of the world's fastest growing economics. It is a 

jrt ' very rich country in terms of natural resources, with oil, diamonds and hydroelectric 

power, and a very fertile agricultural land. However, Angola remains a third world 
country with about one third of its population still depending on subsistence agri- 
culture [30]. In 2009, the life expectancy at birth was 52 years, population median 
age was 17, 45% of the population was under age 15 and only 4% over 60 years, and 
the total fertility rate per woman was 5.6. Angola's total population was 10 million 
(y-\ ' in 1990 and it is approximately double that figure today. According to the Global 

fSJ I Health Observatory Data Repository of the World Health Organization (WHO), in 

CO ■ 2010 the infant mortality rate (probability of dying between birth and age 1) was 

fT^ I 98 per 1,000 live births; and the under- five mortality rate (probability of dying by 

^^ ■ age 5) was 161 per 1,000 live births, which corresponds to 121,000 deaths per year. 

Cn I Tuberculosis (TB) is the leading cause of death among individuals infected with 

the human immunodeficiency virus (HIV) and the rates of co-infection exceed 1,000 

per 100,000 of total population in some countries in sub-Saharan Africa [3]. Be- 

K^ • tween 1990 and 2005, TB incidence rates tripled in African countries with high 

V^ I HIV prevalence [24] . Recent data from Luanda Sanatorium Hospital makes it clear 

C^ i that TB is one of the leading causes of death in Angola [25]. In 2010 alone, there 

were approximately 58,000 new TB cases in Angola. The country established a 

National Program for the Control of Tuberculosis in 1981. However, the civil war 

destroyed 70% of Angola's health facilities, leaving a substantial proportion of the 

population vulnerable to TB. In 2007 the Ministry of Health confirmed that the 

national program was implemented in only 8.6% of health units. Moreover, the 

internal transport and distribution networks for drugs are unreliable and 40% of 

clinics have experienced stockouts in TB drugs [31]. According to WHO's Global 

Tuberculosis Control Report of 2010, the figures of TB in Angola are as follows: 
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the induced mortality rate per 100,000 decreased from 42 in 1990 to 30 in 2009; the 
incidence rate (including HIV) increased from 205 in 1990 to 298 in 2009; the new 
and relapse notification rate increased from 96 in 1990 to 223 in 2009; the number of 
smear-positive cases increased from 3,804 in 2005 to 22,488 in 2009; the number of 
extra-pulmonary TB cases increased from 266 in 1995 to 3,780 in 2010; the relapse 
rate increased from 134 in 1995 to 2,444 in 2010; and the total number of retreated 
individuals increased from 134 in 1995 to 7,776 in 2010. Treatment success rates 
vary between 68% in 2002 and 74% in 2005, and are still within that range today, 
weh below WHO's target of 85% [31, 32]. The multidrug-resistant TB (resistant to, 
at least, isoniazid and rifampicin) was about 2% of all new cases in 2009. 

TB remains a major health problem worldwide. In 2011, the regional director of 
WHO for Africa stated that despite the fact that member states have almost uni- 
versally adopted a strategy to stop tuberculosis, the implementation of this strategy 
has been ineffective. Co-infections TB/HIV and TB multi-resistance to drugs are 
some of the reasons for the difficulty in controlling TB [22] . Most cases of TB are 
caused by the Mycobacterium tuberculosis, which is usually transmitted via airborne 
infection from someone who has active TB. Approximately 10% of infected people 
with Mycobacterium tuberculosis develop active TB disease, that is, approximately 
90% of infected people remain latent. Latent infected TB people are asymptomatic 
and do not transmit TB, but may progress to active TB through endogenous re- 
activation or exogenous reinfection [16, 17]. The anti-TB drugs developed since 
1940 have helped to reduce the mortality rates significantly: in clinical cases, cure 
rates of 90% have been documented and in 2009 Angola's treatment cure rate was 
about 72% [21]. Three types of TB treatment are available: vaccination to pre- 
vent infection; treatment to cure active TB; and treatment of latent TB to prevent 
endogenous reactivation [5]. 

Optimal control is a branch of mathematics that involves finding optimal ways of 
controlling a dynamic system [1, 4, 13, 19]. While the usefulness of optimal control 
theory in epidemiology is nowadays well recognized [2, 9, 10, 11, 14, 15], results 
pertaining to tuberculosis are a rarity [7], and specific studies for the situation of 
Angola nonexistent. Our aim in this paper is thus to use real data from Angola to 
study optimal strategics for the minimization of the number of active TB infectious 
and persistent latent individuals, taking into account the cost of the measures for 
the treatment of these individuals. 

The treatment of active infectious individuals can take different amounts of time 
[8]. Here we consider treatments with a duration of 6 months. With these treat- 
ments one of the barriers to their success is that the patients often do not complete 
them. Since after two months patients no longer have symptoms of the disease and 
feel healed, many of them stop taking the medicines. When the treatment plan is 
not completed, the patients are not cured and reactivation can occur and/or the 
patients may develop resistent TB. The lack of support from family may cause some 
patients to abandon treatment. Statistics from Luanda Sanatorium Hospital show 
that 15 of the 270 TB patients abandoned the hospital during 2010. A possible 
strategy to ensure patients complete the treatment involves supervision and patient 
support. In fact, this is one the measures proposed by WHO's Direct Observation 
Therapy (DOT) [20]. One example of treatment supervision consists of recording 
each dose of anti-TB drugs on the patient treatment card [20] . These measures are, 
however, very expensive, since the patients must stay longer in the hospital or spe- 
cialized people must be paid to supervise patients until they finish their treatment. 
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On the other hand, it is recognized that the treatment of latent TB individuals 
reduces the chances of reactivation, even though it is still unknown how treatment 
influences reinfection [5]. 

In Section 2 we consider a mathematical model for TB from [5] , which considers 
reinfection and post-exposure interventions. We alter this model by adding two 
control functions: the first control ui is associated with preventive measures that 
help active infected patients to complete the treatment; the second control U2 gov- 
erns the fraction of persistent latent individuals that is put under treatment. Our 
aim is stated in precise terms in Section 3: to study how the two control measures 
can reduce the number of infected and persistent latent individuals, taking into 
account the cost associated with their implementation (optimal control problem 
(l)-(3)). Our conclusion is that a successful implementation of these measures can 
prevent the increase of new and reinfected TB cases. In Section 4 we present numer- 
ical simulations for three different strategics and the corresponding optimal control 
solutions are discussed. We end the paper with some conclusions in Section 5. 

2. TB model with controls. We study the mathematical model from [5], where 
reinfection and post-exposure interventions for tuberculosis are considered. To this 
model we add two control functions ui(-) and M2(-) and two real positive model con- 
stants ei and €2- The resulting model is given by the following system of nonlinear 
ordinary differential equations: 

Sit) ^fiN- p{t)Sit) - f,S{t) 
L,{t) = |,/(t) {S{t) + aL2{t) + anR{t)) - (S + r, + ^i) L^{t) 



(1) 



i{t) = (j)5Li{t) + wLaW + i^RR{t) ~ (to + eiui(t) + /i)^(i) 
L2{t) - (1 - (l))5Li{t) - a^I{t)L2{t) -{l^ + e2U2{i) +T2+ ^i)L2{t) 

R{t) = (to + eiui{t))l{t) + TiLi(i) + (t2 + 62^2(^)^2(0 " ^B.j^I{i)R{t) 

^{ujR + fi)R{t). 

The population is divided into five categories (i.e., control system (1) has five state 
variables): 

• susceptible {S); 

• early latent (Li), i.e., individuals recently infected (less than two years) but 
not infectious; 

• infected (/), i.e., individuals who have active TB and are infectious; 

• persistent latent (L2), i.e., individuals who were infected and remain latent; 

• and recovered (R), i.e., individuals who were previously infected and have 
been treated. 

The control ui represents the effort in preventing the failure of treatment in active 
TB infectious individuals /, e.g., supervising the patients, helping them to take 
the TB medications regularly and to complete the TB treatment. The control U2 
governs the fraction of persistent latent individuals L2 that is put under treatment. 
The parameters e^ g (0,1), i = 1,2, measure the effectiveness of the controls Ui, 
i = 1,2, respectively, i.e., these parameters measure the efficacy of treatment inter- 
ventions for active and persistent latent TB individuals, respectively. We assume 
that the total population, N, with N = S{t) + Li{t) + I{t) + L2(t)+R{t), is constant 
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in time. In other words, we assume that the birth and death rates are equaL By 
virtue of this assumption, we can reduce the control system (1) from five to four 
state variables. However, in this paper we prefer to keep the TB model in form (1), 
using relation S (t) + Li{t) + 1 (t) + L2{t) + R{t) ^ N a.s a test to check the numerical 
results. 

Following who's Global Health Observatory for Angola in 2009 [23], we consider 
IJ, = l/52yr~^, which corresponds to a life expectancy at birth of 52 years. This 
and other values of the parameters are given in Table 1. Moreover, following [5], we 
assume that there are no disease-related deaths, and that at birth all individuals are 
equally susceptible and become different as they experience infection and respective 
therapy. The proportion of the population in each category changes according to 
the system (1). 

Two different values for the initial conditions S{0), Li{0), /(O), ^2(0) and R{0) 
are considered (sec Table 2): 

• Case 1, initial values are borrowed from [7]; 

• Case 2, initial values are based on the data from Luanda's Sanatorium Hos- 
pital during September 2010 [25]. 

The values of the rates S, (f), uj, lvj^, a and tq arc taken from [5] and the references 
cited therein (see Table 1 for the values of all the parameters). The parameter S 
denotes the rate at which individuals leave the Li category; (/) is the proportion of 
individuals going to category /; uj is the rate of endogenous reactivation for persis- 
tent latent infections (untreated latent infections) ; and ojr is the rate of endogenous 
reactivation for treated individuals (for those who have undergone a therapeutic in- 
tervention). The parameter cr is a factor that measures the reduction in the risk of 
infection, as a result of acquired immunity to a previous infection, for persistent la- 
tent individuals, i.e., this factor affects the rate of exogenous reinfection of untreated 
individuals; while an represents the same parameter factor but for treated patients. 
In our simulations we consider the case where the susceptibility to reinfection of 
treated individuals equals that of latents: aR = a. The parameter tq is the rate of 
recovery under treatment of active TB (assuming an average infectiousness duration 
of six months). The parameters ti and T2 apply to latent individuals Li and L2, 
respectively, and are the rates at which chemotherapy or a post-exposure vacine 
is applied. In [5], different values for these rates are considered: the case where 
no treatment of latent infections occur (ti ~ T2 =0); the case where there is an 
immediate treatment of persistent latent infections (t2 — >■ cx)); and the case where 
there is a moderate treatment of persistent latent infections (t2 S [0.1,1]). The 
first and second cases are not interesting from the optimal control point of view. In 
our paper we consider, without loss of generality, that the rate of recovery of early 
latent individuals under post-exposure interventions is equal to the rate of recovery 
under treatment of active TB, ti = 2 yr^^, and greater than the rate of recovery of 
persistent latent individuals under post-exposure interventions, T2 = lyr~^. It is 
assumed that the rate of infection for susceptible individuals is proportional to the 
number of infectious individuals and the constant of proportionality is /3, which is 
the transmission coefficient. 

From the cpidcmiologistic point of view, an important quantity is given by the 
basic reproduction number Rq. This is the expected number of secondary cases 
produced in a complete susceptible population, by a typical infected individual 
during his/her entire period of infectiousness [15]. In the absence of controls (ui = 
U2 = 0), it is known (see [5]) that the basic reproduction number for system (1) is 
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proportional to the transmission coefficient /3: 



Ro = P 



^i{llib. + To + fi){5 + fi)(uj + ^) 



The endemic threshold is given at i?o = 1 and indicates the minimal transmission 
potential that sustains endemic disease, i.e., when Rq < 1 the disease will die out 
and for Rq > 1 the disease may become endemic. We can prove that variations of 
the parameter /3 are associated with different values for the optimal control ui but 
do not interfere significantly with the optimal control U2- Moreover, as expected, 
as (3 decreases the fraction of latent and infected individuals also decreases. Here 
we consider (3 — 100, which corresponds to the endemic case {Rq ~ 2.96 > 1). 
The measures of control efficacy are fixed to the values e^ = 1/2, i = 1,2. Because 
the risk of developing disease after infection is much higher in the first five years 
following infection, and declines exponentially after that [18, Sec. 5.5], we take the 
total simulation time duration as T = 5 years. 



3. The optimal control problem. Consider the state system (1) of ordinary 
differential equations in M.^ with the set of admissible control functions given by 

n = {(ui(-),u2{-)) e (L°°(o,r))2 1 o < Ui{t),u2{t) < i, vt e [o,r] }. 

The objective functional is given by 

J{ui{-),U2{-))^ f [l{t)+L2{t) + ^Wiul{t) + ^W2ul{t)]dt, (2) 

Jo 

where the constants Wi and W2 are a measure of the relative cost of the interventions 
associated with the controls ui and U2, respectively. In our numerical simulations 
we consider the weight Wi associated with the control ui to be greater than the 
weight W2 associated with U2: W\ = 500 and W2 = 50. The reason for this is that 
the cost associated with ui includes the cost of holding active infected patients / 
in the hospital or paying professionals to supervise them, ensuring that they finish 
their treatment, which is costly to implement. In Section 4.4 we compare this situa- 
tion to the one when the weights take the same values: Wi = W2 = 50. We consider 
the optimal control problem of determining {S*{-),Ll{-),I*{-),L2{-),R*{-)), associ- 
ated with an admissible control pair (ui(-); "2(')) G f^ on the time interval [0,T], 
satisfying (1), the initial conditions S'(O), Li(0), J(0), L2(0) and R{0) (see Table 2), 
and minimizing the cost functional (2), i.e., 

JK(-),^(-))-nnnJK(.),U2(-)). (3) 

Theorem 3.1. The problem (1), (3) with fixed initial conditions S{0), Li{0), 
/(O), ^2(0) and R(0), and fixed final time T, admits a unique optimal solution 
{S* {■) ^ L\{-) , I* {■) ^ L2{-) T R* {■)) associated with an optimal control pair {u\{-) , U2{-)) 
on [0,T]. Moreover, there exist adjoint functions, A|(-), A2(-), A3(-), A4(-) and 
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A5(-), such that 



Kit) = Kit){j^i*it) 



m)j^i*it) 



X*{t) = X*{t)iS + n+ aO - A,;(i)0<5 - A:(i)(l - 0)<5 - A^(i)ri 



13 



li 



Xm = -1 + XUt)^S*{t) - xmT?iS*it) + aLl{t) + anR*{t)) 



N 



N 



+A*(t)(ro + 6i<(i) + m) + A:(i)^^i2(t) - Kit) [t^ + e,ul{i) - aR^R*{i) 

Xlit) = -1 - A*(i)-^r (t)a - A*(t)w + A^W ff^^^* W + ^ + ^2u;{t) +r2+fi 
-Xlit)iT2 + e2u;{t)) 

KW = -A;(t)cTfl|rr W - xm^B + Kit) (oR^i*it) + uJR + ^l 



with transversality conditions X*{T) — 0, i — I, ... ,5. Furthermore, 

eiriK-K) 



<it) 



min < max < 0. 
min < max < 0, 



Wi 
e^LlJXl - Xi: 

Wo 



,1 



(4) 



(5) 



Proof. The existence of optimal controls (u*(-), M2(')) ^"^^ associated optimal solu- 
tion [S* , L\,I* ^L\,R*) comes from the convexity of the integrand of the cost func- 
tion (2) with respect to the controls (wi, U2) and the Lipschitz property of the state 
system with respect to state variables {S,Li,I,L2,R) (see, e.g., [1, 4] for existence 
results of optimal solutions). According to the Pontryagin Maximum Principle [13], 
if (w^(-), W2(')) G ^ is optimal for the problem (1), (3) with the initial conditions 
given in Table 2 and fixed final time T, then there exists a nontrivial absolutely 



continuous mapping A : [0,r] — > 
the adjoint vector, such that 



A(t) = (Ai(i),A2(i),A3(i),A4(t),A5(i)), called 



and 



s 


dH 
dXi ' 


ii 


dH 

0X2' 




dH 


< / 


8X3 ' 




dH 


L2 


dXi' 


R 


dH 

dX5 


'Ai = 


dH 

dS ' 


A2 = 


dH 

dLi ' 


A3 = 


OH 

ai ' 


A4 = 


dH 
dL2 ' 


A5 = 


dH 
dR ' 



(6) 
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where the Hamiltonian H is defined by 
H = H{S, Li,I, L2, R, A, Ml, -U2) 

= / + L2 + — Ul + — "2 

+ Ai LlN - ^/5 - ^iS 

+ A2 f -^/ (5 + (TL2 + ariR) -{5 + ti+ fi)Li 
+ A3 {(f>SLi + UJL2 + ujrR - (to + (-lUi + ii)I) 
+ A4 ( (1 - (l))6Li - a—IL2 - (w + e2U2 + T2 + n)L2 

+ A5 ( (tq + eiMi)/ + TiLi + (t2 + £2^2) L2 - (Tr—IR - (yjR + [i)R 

and the minimization condition 

H{S\t),L\(t),r{t),L*2{t),R*{t),\*(t),u\(t),u*2{t)) 

= min H{S*{t),Ll{t),r{t),L'2it),R*it),X*{t),ui,U2) (7) 

0<«1,M2<1 

holds almost everywhere on [0,r]. Moreover, the transversality conditions Xi{T) = 
0, i = 1, . . . , 5, hold. System (4) is derived from (6), and the optimal controls (5) 
come from the minimization condition (7). For small final time T, the optimal 
control pair given by (5) is unique due to the boundedness of the state and adjoint 
functions and the Lipschitz property of systems (1) and (4) (see, e.g., [7]). D 

4. Numerical results and discussion. In this section we present different op- 
timal control strategies for the TB model (1) under the parameter values given in 
Table 1. Three control strategies are explored: 

• Strategy 1, which implements measures that prevent the failure of treatment 
in active TB individuals / (control ui alone); 

• Strategy 2, which considers treatment of persistent latent individuals L2 
with anti-TB drugs (control 1*2 alone); 

• Strategy 3, which implements measures for preventing treatment failure in 
active TB individuals and variation of the fraction of persistent latent indi- 
viduals that are put under treatment (controls ui and U2). 

Different approaches were used to obtain and confirm the numerical results. One 
approach consisted in using IPOPT [26] and the algebraic modeling language AMPL 
[27]. A second approach was to use the PROPT Matlab Optimal Control Software 
[28]. The results coincide with the ones obtained by an iterative method that 
consists in solving the system of ten ODEs given by (1) and (4). For that, first 
one solves system (1) with a guess for the controls over the time interval [0,T] 
using a forward fourth-order Runge-Kutta scheme and the transversality conditions 
Xi{T) = 0, i — 1,...,5. Then, system (4) is solved by a backward fourth-order 
Runge-Kutta scheme using the current iteration solution of (1). The controls are 
updated by using a convex combination of the previous controls and the values from 
(5). The iteration is stopped if the values of unknowns at the previous iteration are 
very close to the ones at the present iteration. For more details see, e.g., [7]. 
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Symbol 


Description 


Value 


P 


Transmission coefficient 


100 


M 


Death and birth rate 


1/52 j/r'^ 


5 


Rate at which individuals leave Li 


12 yr"^ 


^ 


Proportion of individuals going to / 


0.05 


LU 


Endogenous reactivation rate for persistent latent infections 


0.0002 yr-i 


UJR 


Endogenous reactivation rate for treated individuals 


0.00002 yr^^ 


a 


Factor reducing the risk of infection as a result of acquired 






immunity to a previous infection for L2 


0.25 


OR 


Rate of exogenous reinfection of treated patients 


0.25 


TO 


Rate of recovery under treatment of active TB 


2yr-^ 


Tl 


Rate of recovery under treatment of latent individuals Li 


2yr-i 


T2 


Rate of recovery under treatment of latent individuals L2 


lyr-^ 


N 


Total population 


30, 000 


T 


Total simulation duration 


5 yr 


ei 


Efficacy of treatment of active TB I 


0.5 


£2 


Efficacy of treatment of latent TB L2 


0.5 


Wi 


Weight constant of control ui 


500 


W2 


Weight constant of control U2 


50 



Table 1. Parameter values for the optimal control problem (l)-(3) 



States 
(initial value) 


Case 1 
(fraction) 


Case 1 

(%) 


Case 2 
(fraction) 


Case 2 

(%) 


5(0) 


(76/120)A/- 


~ 63.3% 


(16/50)iV 


32% 


Li(0) 


(36/120)Af 


30% 


(28/50)iV 


56% 


7(0) 


(5/120)N 


~ 4.2% 


(3/50)Af 


6% 


L2{0) 


(2/120)iV 


~ 1.7% 


{2/50)N 


4% 


R{0) 


(l/120)iV 


~ 0.8% 


{1/50)N 


2% 



Table 2. Initial values for the state variables of (1) 



To begin, we study the influence of the initial values of the state variables, as 
given in Table 2, on the optimal solution associated with Strategy 3. 

4.1. DifTerent initial values (Case 1 and Case 2) for the state variables. 

Our aim is to compare the optimal solution associated with Strategy 3 (when both 
controls ui and U2 are considered) under different initial values for the state vari- 
ables: in the first case we consider the initial values of [7] (Case 1) and in the second 
case we consider the initial values given in an official communication from Luanda 
Sanatorium Hospital [25] (Case 2). In Figure 1 we observe that the differences 
between the initial values considered in Cases 1 and 2 (Table 2) do not influence 
the optimal controls significantly. In Case 1, the initial percentage of susceptible 
individuals 5(0) is approximately twice the respective percentage in Case 2. How- 
ever, with the optimal controls (wj, Uj) illustrated in Figure 1, after two years this 
difference is substantially reduced, as we can observe in Figure 2. With respect to 
the initial values of early latent individuals Li{0), the percentage in Case 1 is 30% 
and in Case 2 it is 56%. Here, in both cases, the percentage of early latent indi- 
viduals is reduced to approximately 10% after approximately 140 days. We observe 
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(a) Optimal Control u^ 



(b) Optimal Control tij 



Figure 1. Optimal control functions in Case 1 (solid line) and 
Case 2 (dashed line), in a time scale in units of years, for Strategy 3 
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0.3 
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0.2 


! 


0.1 


/ Case 1 ■ 




/ Case 2 







(C) Ll/N 



(d) R'/N 



Figure 2. Optimal states in Case 1 (solid line) and Case 2 (dashed 
line), in a time scale in units of years, for Strategy 3 



that even if in Case 2 the initial fraction of early latent individuals is much higher 
than in Case 1, the optimal controls ensure a big reduction in the fraction of early 
latent individuals and, after no more than 5 months, the difference between the 
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two cases no longer exists. Given these results, from now on we consider, without 
loss of generality. Case 1 for the initial values of the state variables. Our numerical 
simulations were checked using the relation S{t) + Li{t) + I{t) + L2{t) + R{t) ~ N . 
Figure 3 shows that S{t) + Li{t) + I{t) + L2{t) + R{t) is constant along time for 
Strategy 3, as expected. 




Figure 3. Plot of {S{t) + Li{t) + I{t) + L2{t) + R{t)) /N versus 
t in years for Strategy 3 



4.2. Strategy 1 versus Strategy 3. If we choose to apply only the control ui 
(Strategy 1) and compare this situation to the case where both controls ui and 
U2 are applied (Strategy 3), we observe that the optimal control u\ stays at the 
upper bound for almost the same duration in both situations (sec Figure 4). When 
only the control ui is considered, the fraction of persistent latent individuals L2 is 
greater than the respective fraction associated with Strategy 3, for the entire five 
years (see Figure 5 A and E) . This implies a higher value of the cost functional (2) 
associated with Strategy 1 when compared to the cost associated with Strategy 3. 
When one compares the change in the fraction of susceptible (5"), early latent (Li) 
and infected (/) individuals, no difference is observed between Strategies 1 and 3 
(see Figure 5 B, C and D). 
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(a) Control mJ 

Figure 4. The optimal control u\ for Strategy 1 (dashed line) and 
Strategy 3 (solid line) in a time scale in units of years 
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(A) (7* + L^)/7V 
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(B) S*/N 
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(E) Ll/N 



(f) R'/N 



Figure 5. Optimal states for Strategy 1 (dashed line) and Strat- 
egy 3 (solid line) in a time scale in units of years 



4.3. Strategy 2 versus Strategy 3. When we compare Strategy 2 (only control 
U2) with Strategy 3 (both controls ui and M2), wc observe that the best choice is 
definitely to use Strategy 3. Indeed, with Strategy 3, there is a lower fraction of 
early and persistent latent individuals as well as infected individuals (see Figure 6 
A, C, D and E). The fraction of susceptible and recovered individuals is higher when 
two controls are applied at the same time (see Figure 6 B and F). Similarly with 
Strategy 1, the value of the cost functional (2) associated with Strategy 3 is lower 
than that of Strategy 2. 
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Figure 6. Optimal states for Strategy 2 (dashed line) and Strat- 
egy 3 (solid line) in a time scale in units of years 



4.4. Changing the relative cost of interventions. So far all the simulations 
were done assuming that the weight Wi associated with the control wi is greater 
than the weight W2 associated with U2: Wi = 500 and W2 = 50 (see Table 1). We 
now consider the case when the weights are the same: Wi = W2 = 50. In Figure 7 
we observe that a higher weight associated with the control ui implies that the 
optimal control u* stays at the upper bound for a smaller duration. Surprisingly, 
this change in the control does not result in a change of the behavior of the state 
variables (see Figure 8). The results with Wi = W2 = 50 for Strategy 1 are similar 



OPTIMAL CONTROL STRATEGIES FOR TUBERCULOSIS 



13 




(a) Optimal control u^ 

Figure 7. Optimal control u* for Strategy 3 in a time scale in 
units of years: Wi = 500 and W2 = 50 (solid line) versus Wi = 
W2 = 50 (dashed line) 



to the ones reported in Section 4.2 (including the number of days that the optimal 
control u* stays at the upper bound). 

5. Conclusion. Tuberculosis (TB) is an important health issue all over the world, 
particularly in many African countries. In this article we focused our attention in 
Angola where, since 2009, it is recognized by authorities that TB and HIV are the 
two most serious public health problems. Efforts to control TB in Angola have 
been underway since 1981, but results are still not satisfactory, and minimizing the 
effects caused by TB is an important challenge. Implemented measures to control 
these problems are having a positive effect, but the reduction rate is very slow [29]. 
We introduced two control functions ui and U2 to an existing mathematical 
model for TB developed in [5]. These controls are associated with measures that 
help to reduce the number of active infected and persistent latent TB individuals: 
the control wi represents the effort that prevents the failure of treatment in active 
TB infectious individuals /, e.g., supervising the patients, helping them to take the 
TB medications regularly and to complete the TB treatment; while the control U2 
governs the fraction of persistent latent individuals L2 under treatment with anti- 
TB drugs. An optimal control problem was formulated and solved theoretically 
using the Pontryagin maximum principle. The solution to the problem was then 
illustrated by numerical simulations using available data from Angola. From the 
numerical results, one may conclude that considering only control ui (Strategy 1) 
or only control U2 (Strategy 2) does not lead to the best results in terms of the 
number of active infected and persistent latent individuals. A combined strategy 
(Strategy 3) that involves both controls is preferable. Figure 9 shows the simulation 
of the system for the uncontrolled case, when ui = and U2 = 0, compared with the 
strategies here proposed. Results justify the need for intervention in tuberculosis 
treatment. The values of Table 3 put in evidence that the uncontrolled situation 
is the worst: does not only result in more infected and persistent latent individ- 
uals (1,584 individuals at the end of five years versus 334 individuals in the case 
Strategy 3 is applied) but also gives a higher value to the cost functional (2). 
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Strategy 1 


Strategy 2 


Strategy 3 


Without controls 


/(5) + i2(5) 


-751 


-831 


-334 


- 1, 584 


Cost: functional (2) 


~ 32, 511.2 


-28,585.9 


-24,133.1 


- 37, 760.3 



Table 3. Comparison of Strategies 1, 2, 3 with the uncontrolled 
case, i.e., when ui = M2 = 
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